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ABSTRACT 

We present high resolution submilhmeter interferometric imaging of two of the brightest high- 
redshift submilhmeter galaxies known: GN20 and AzTECl at 0.8 and 0.3 arcsec resolution respec- 
tively. Our data - the highest resolution submilhmeter imaging of high redshift sources accomplished 
to date - was collected in three different array configurations: compact, extended, and very ex- 
tended. We derive angular sizes of 0.6 and 1.0 arcsec for GN20 and 0.3 and 0.4 arcsec for AzTECl 
from modeling their visibility functions as a Gaussian and elliptical disk respectively. Because both 
sources are B-band dropouts, they likely lie within a relatively narrow redshift window around z ~ 4, 
which indicates their angular extent corresponds to physical scales of 4-8 and 1.5-3 kpc respectively 
for the starburst region. By way of a series of simple assumptions, wc find preliminary evidence that 
these hyperluminous starbursts - with star formation rates > 1000 Mq yr~^ - are radiating at or close 
to their Eddington limit. Should future high resolution observations indicate that these two objects 
are typical of a population of high redshift Eddington-limited starbursts, this could have important 
consequences for models of star formation and feedback in extreme environments. 
Subject headings: cosmology: observations - galaxies: evolution - galaxies: high-redshift - galaxies: 
starburst - galaxies: submilhmeter - galaxies: formation 



1. INTRODUCTION 



Wide area surveys at n iillim eter (e.g., iGreve et al.l 
12004 iBertoldi et all 120071: IScott ct al. 2008) an d sub- 
millimeter (e.g., Small et al.l 11997.: iHughcs et al.l[l998l : 



iBarger et al]|1998l : lPope et al.ll2006l : ICoppin et al.]|200e 

wavelengths have revealed a large population of ultra- 
and hyperluminous infrared galaxies - ULIRGs and 
HyLIRGs - at high redshift (median z ~ 2 for a radio- 
selected sample [chapman et al.l[2005l ). Since their initial 
discovery it has become clear that these "submilhmeter 
galaxies" (SMGs) are likely massive, g as-rich merging 
svstems (iFraver et al. 1998, 1999 ; Chap man et~al]|2003t 
IGreve et al.l l2005t iTacconi et al.l [200 6. 2008) tha t rep- 
resent massive g alaxies in formation (,Scott et al.ll2002l : 
iBlain et al.l 12004 ). Extremely luminous infrared objects 
take on in creasing cosmologica l importance at z ^ 1 
([Sanders fc Mhabcl 1996; Lc Floc'h et all [200^), and 
may dominate cosmic star formation for up to the first 
half of the lifetime of the universe ([Blain et al.l Il999l . 
[200l . 

Despite significant progress over the past decade, a 
more complete understanding of SMGs has been ham- 
pered in part by the relatively poor resolution of submil- 
hmeter cameras [^ 10 — 18 arcsec FWHM). In particular, 
the size scale of the starburst region - traced by the rest 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138 

^ jyounger@cfa.harvard.edu 

^ Institute of Astronomy, The University of Tokyo, 2-21—1 
Osawa Mitaka, Tokyo, 181-0015, Japan 

* Current Address: National Radio Observatory, NAOJ, Mi- 
namimaki, Minamisaku, Nagano, 384-1305, Japan 

5 Joint ALMA Office, El Golf 40, Las Condes, Santiago 7550108, 
Chile 

® Astronomy Department, University of Massachusetts, 
Amherst, MA, 01003 



frame far infrared (IR) - potentially provides important 
insights into the nature of the engine driving the tremen- 
dous luminosity of these systems. If they are scaled up 
versions of local ULIRGs, we would expect far-IR emis- 
sion on scales of ~ 5 — 10 kpc (e.g.. lDownes fc SolomonI 
[19981: llono et al.ll2007l see also lono et al. 2008, in prepa- 
ration). Hydrodynam ic modeling of merger dri ven nu- 
clear starbursts (e.g., [Mihos fc Hernquiff '1994") of the 
kind thought to drive many SMGs ((Chapman et al.l 
[20031: [Tacconi et all [20061 [2008[ ) can be somewhat more 
compact, a result that could have important physical 
consequences: Eddington argume nts suggest a mini- 
mum s ize scale for such regions ([Murrav et al.[ [2005[ : 
[Thomp son et al. 2005). Unfortunately, at typical SMG 
redshifts, all of these size scales are far smaller than the 
typical angular resolution of submillimeter cameras on 
single-dish instruments. 

The first breakthrough came with deep radio contin- 
uum surveys, which leveraged the local far-IR/radio cor- 
relation (iCondonI I1992D in combinat ion with statistical 
arguments ([Ivison et al.l [2002L [20071 ) to associate faint 
radio counterparts within the submillimeter beam with 
SMG s. Higher resoluti o n radio imaging of the se sources 
(Cha pman et~al] [2004 [Biggs fc IvisonI [20081) found a 
range of source structures on physical scales of '^ 1—8 kpc 
with a median of 5 kpc. While promising, this technique 
assumes a spatially resolved far-IR/radio correlation 
which is not particul a rly well t inders t ood locally (e.g., 
Hippelein et~all [2003[ : [Murphvl [200l [Tabatabaei et all 



2003) ■ As a consequence, these results are not straight- 



forward to interpret. 

This motivates high-resolution imaging of the rest 
frame far-IR directly, via submillimeter interferometry. 
The vast majority of previous work was done at res- 
olutions of '--^ 1 — 2 arcsec, and found that the far- 



TABLE 1 
Track Details 



Target 


Configuration u — v Coverage 
[kA] 


Beam Size 
[arcsec] 


Date 
[dd.mm.yy] 


< T225GHZ > 


Obs. Time" 

[hrs] 


Reference'' 


GN20 


COM 
EXT 


15-75 
40-200 


2.99 X 2.26 
0.81 X 0.75 


20.02.05, 05.03.05 
10.02.08, 11.02.08 


0.04, 0.06 
0.04, 0.04 


10.4 
5.3 


106 
This work 


AzTECl 


COM 
EXT 
VEX 


20-75 
50-250 
60-550 


2.69 X 2.19 
0.86 X 0.55 
0.25 X 0.35 


17.01.07 
16.01.08 
05.04.08 


0.05 
0.03 
0.03 


5.6 
4.0 
3.4 


Y07 
This work 
This work 



Total on— source integration time in that configuration.'' 106: ,Iono et al.l II2006I ): Y07: I Younger et al.l II2007I ) 
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Fig. 1. — The u — v coverage for our high resolution interferometric imaging of GN20 (left) and AzTECl (right). Included are all tracks 
in three different SMA configurations: compact (COM: black), extended (EXT: blue), and very extended (VEX: red). For further details, 
including weather conditions and on-source integration times, see Table [Tl 

IR continuum in SMGs orig inates at physical scales 
of ^4-8 kpc |Neri et al.i 12 003: Grove et al. 120051: 
Tacconi et ap l200a iWang et al.. ,2007.: .Younger et all 
20071 120081 : iDannerbaue r et al.1 I2008D. More re cently. 



very high resolution CO imaging bvlTacc oni et al.1 ()2008D 
showed that gas motions in typical SMGs are disordered 
on scales of ~ 1 — 2 kpc, suggesting that they are ongoing 
major mergers. 

In this paper, we present high resolution (beam size 
^ 1 arcsec) 890 /im continuum imaging of two of 
the bright e st SM Gs known - GN20 (Pope ct al. 2006 



lono et aT] 



200(| and AzTECl (lYouneer et all l20()7l : 



Scott et al]|2008l) - with the Submillimeter Array (SMA 
Ho et al.ll2004[ ). By targeting t he brightest - and there- 



fore likely the most luminous ("Blai n fc Longaiiill993( ) 
objects, we constrain the physical scale of the far-IR 
in extreme conditions. In addition, since these objects 
are thou ght to lie at higher re dshift than radio-selected 
samples () Younger et al.|[2007f) . they offer an intriguing 



probe of the nature of star formation at earlier epochs. 
Throughout this work we assume a concordance cosmol- 
ogy with rim = 0.3, fl\ — 0.7, and h = 0.7. 

2. OBSERVATIONS AND DATA REDUCTION 

The two targets were GN20 - the brightest 850 /im 
source in the S ubmillimeter Comrn on-User Bolometric 
Array (SCUBA: [Holland et anfl999l) survey of the Hub - 
ble Deep Field North (HDFN: see iPope et all l2006f ) . 
and AzTECl - the brightest 1.1mm source in the 
AzTEC (Wilson et al.l[2008:') survey of the COSMOS field 
(jScott et al.ll2008f) . Both were previously detected as sin- 



gle point-sources wi th flux densities of J sgoum = 22.9 ± 
2.8 (llono et al.ll2006D and 15.6 ±1.1 mJy (jYounger et all 
|2007| ) respectively with the SMA in compact configu- 
ration (COM). We have re-observed both these targets 
with the SMA in extended configuration (EXT), which 
provides a ~ 3x improvement in angular resolution over 
COM, using the same pointing center as the COM tracks. 
The EXT tracks - two for GN20, one for AzTECl - were 
taken in excellent weather in January and February 2008. 
Since AzTECl was unresolved in the EXT track (see §[3] 
and Figure [3|), we re-observed it in very extended config- 
uration (VEX) in April 2008, which provided a further 
'^ 3x improvement in angular resolution. For details on 
the tracks, configurations {u — v range, beam size, etc.), 
and observing conditions, see Table [TJ 

The receiver was tuned to 345 GHz in the USB, and 
averaged with the LSB for an effective bandwidth of 4 
GHz centered at 340 GHz. For GN20, passband calibra- 
tion was done using 3C273 and 1921-293, and primary 
flux calibration was done using Titan. The target was ob- 
served on a 10 minute cycle - 5 minutes on source, 5 min- 
utes on calibrators - with two primary gain calibrators: 
1048-f717 (~ 0.3 Jy; 14 degrees away) and 1153-^495 
(~ 0.3 Jy; 14 degrees away). For AzTECl, passband cal- 
ibration was done using 3C111 and 3C273, and primary 
flux calibration was done using Ceres. As with GN20, 
the target was observed on a 10 minute cycle with two 
primary gain calibrators: 1058-1-015 (~ 2 Jy; 15 degrees 
away) and 0854-1-201 (~ 2 Jy; 24 degrees away). Because 
Ceres is known to be variable at the ^^ 20 — 30% level due 
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Fig. 2. — Stamp images, 3 arcsec on a side, of GN20 (top) and AzTECl (bottom): (from left to right) SMA 890 /im dirty map using 
EXT configuration data, B-band imaging (HST/ACS for GN20, ground based Subaru for AzTECl), i-band HST/ACS imaging, and VLA 
20cm imaging data. The contours overlaid on all frames are from the SMA EXT configuration dirty map, in intervals of 3,5,7... times 
the r.m.s. noise. For reference, the grey dashed line is 1 arcsec in length and the beam is indicated for each with a hashed grey ellipse: 
0.81 X 0.75 and 0.86 x 0.55 arcsec FWHM for GN20 and AzTECl respectively. 



to rotation (lAltenhoff et all 119941 : iRedman et"aLl 119981 : 
iBarrera-Pineda et al.l |2005[) . we confirm this fiux scale 
by checking that the flux density for 0854+201 derived 
from this track (F340GHZ = 2.37 Jy) is consistent to that 
measured one day earher (F340GHZ = 2.29 ± 0.12 Jy) for 
which Titan was the primary flux cahbrator. 

In addition to the two primary targets, we observed 
a nearby test quasar once every 60 minutes through- 
out the track to empiricaUy verify the phase transfer 
and inferred source structure, and estimate the sys- 
tematic positional uncertainty. The test quasars for 
GN20 and AzTECl were Jl'302+578 (- 0.1 Jy; 5.5 
degrees away) and J1008+063 (~ 0.2 Jy; 5 degrees 
away) respectively. Both are included in both the JVAS 
(iPatnaik et al.lll992 ':'Brownc et al."1998') and VLB A Cal- 
ibrator (|Ma et all 11998: Bcasley ct al. 2002) surveys of 
compact, flat-spectrum radio sources, and have absolute 
positions known to better than 20 mas. 

For the VEX track on AzTECl, time dependent gain 
calibration dervied from 1058-1-015 left clear, slow, resid- 
ual phase variations on 0854-f 201 (and J1008+063) due 
to uncertainties in the baseline parameters or other lim- 
itations of the SMA interferometer model. To improve 
the phase transfer and prevent decor relation, an addi- 
tional gain calibration was performed using J1008-I-063, 
just 4 degrees away from AzTECl in declination, since 
conditions were good enough to yield sufficient signal-to- 
noise (> lOcr) in each of the hourly scans on this source. 
This additional step minimized the phase errors owing to 
baseline effects in the calibrated visibilities. Remaining 
phase fluctuations dominated by the atmospheric effects 
on short timescales left an effective seeing sizescalc of 
- 0.08 arcsec in J1008-I-063 (see § [3] and Figure g] for 
further discussion). 

We also make use of extensive multiwavelength data in 
both fields. For the HDFN, thi s includes HST/ACS B - 
, V-, i-, and z-band optical (Giavalisco et al.ll2004af). 
IRAC 3.6-8.0 iim and MIPS 24 ;nm (jDickinson et all 



I2003f l. and VLA 20cm (iBiggs &: IvisonI l2006tl imag- 
ing d ata. For the COSMOS field (see IScoville et all 
I2OO7I for an o verview), this includes Subaru ground 
based optical (iTaniguchi et all l2007f l. HST/ACS i- 
band (jKoekem oer et al.' '20(0 IRA C 3.6-8 fim and 
MIPS 24 nm (SandCTS et aTr i2007|), and VLA 20cm 
(jSchinnerer et al.ll2007t ) imaging. 

3. RESULTS 

Both targets were detected at high significance by the 
SMA in EXT configuration with a ~ 0.75 arcsec beam. 
The maps, along with overlays on multiwavelength imag- 
ing data, are presented in Figure [21 Source structure de- 
rived from the calibrated visibilities - which show fiux 
density as a function of decreasing angular scale - and 
the empirical verification of phase transfer are summa- 
rized in Figures [3] and m and in Table [2 

GN20 shows evidence of being partially resolved by the 
SMA in EXT configuration, with a characteristic angu- 
lar scale of ^ 0.5 — 1.2 arcsec (see Figure [3]and Table[2]) 
as inferred from modeling its visibility as both a Gaus- 
sian and elliptical disk. Its submillimeter position is co- 
incident with a bright IR AC 3.6-8 Mm and faint MIPS 
24 ^m (-F24/jm ^ 70 ^Jv: iPope et al.l [20061) source, and 
roughly so with a radio source (i^20cm = 57 ± 10 /iJy; 
^ 0.5 arcsec away). The SMA map is also consistent 
with the radio morphology, which shows some evidence 
of being resolved along its major axis w ith a beam size of 
1.5 X 1.5 arcsec (jBiggs fc Ivison|[2006[) . High-resolution 
(0.08 arcsec PSF) ACS imaging shows that the submil- 
limeter detection is not coincident with the nearby opti- 
cal "smudge" . However, since this source is a B-band 
dropout, which suggests a redshift range c onsistent with 
the observed radio-to-sub millimeter (see ICarilli fc YunI 
119991: lYun fc Carillil 12(3021) and 24um-to-s ubmillimeter 
(see lWang et al.ll2007l : lYounger et al.ll2007f l fiux density 
ratios of GN20 (see §[4|), it is plausible that this object is 
physically associated with GN20 and represents a region 
of lower dust opacity. 



TABLE 2 
Positions and Source Structure 



Name 



Config." 



Model 



[J2000] 



5 
[J2000] 



farcsecl 



farcsecl 



-FsgoMm 

[mjy] 



maj 

[arcsec] 



min 

farcsecl 



[deg] 



GN20 



C 

E 

E 
C+E 

C+E 



Point 
Gaussian 

Disk 
Gaussian 

Disk 



12:37:11.920 
12:37:11.898 
12:37:11.897 
12:37:11.903 
12:37:11.901 



+62:22:12.17 
+62:22:12.14 
+62:22:12.16 
+62:22:12.16 
+62:22:12.17 



0.10"= 
0.06 
0.06 
0.06 
0.05 



0.10= 
0.09 
0.09 
0.09 
0.06 



22.9 + 2.8 
26.9 + 5.1 
27.3 + 4.5 
23.9 + 2.6 
24.2 + 2.5 



0.8 + 0.2 
1.4 + 0.3 
0.8 + 0.2 
1.2 + 0.3 



0.5 + 0.2 
0.7 + 0.3 
0.3 + 0.3 
0.8 + 0.3 



15 
25 
35 
-25 



AzTECl 



E 
C+E 

yf 
yf 



Point 
Point 
Point 
Gaussian 
Disk 



15.6 + 1.1 
13.8 + 2.3 
15.1 + 1.1 
15.1 + 1.1 
15.1 + 1.1 



09:59:42.859 
09:59:42.863 
09:59:42.863 
09:59:42.863 
09:59:42.863 



+02:29:38.21 
+02:29:38.19 
+02:29:38.19 
+02:29:38.20 
+02:29:38.20 



0.11 
0.07 
0.05 
0.04 
0.04 



0.20 
0.07 
0.06 
0.06 
0.06 



'0.3 

-0.4 



0.2 
0.3 



40 
25 



'^ Data restricted to tliis combination of configurations: COM (C), EXT (E), and VEX (V). See Table [T] and Figure [T] for 
details.'' Combined statistical and systematic uncertainty, where the systematic uncertainty is estimated from the position of 
the test quasar.'^ f'maj a-nd ^min represent the FWHM or diameter of the major and minor axes for the Gaussian and elliptical disk 
models respectively.'* Position angle." Since there was no test quasar available for this track, these are just the statistical positional 
uncertainties from llono et al. (2006).' Because there was some residual phase error due to atmospheric seeing (see §[2]for details), 
we chose not to combine the V track results with C and E, which did not show these errors. Furthermore, we fixed the total flux 
to that inferred from the CE track. The ~ indicates that these size measurements are large compared with the effective seeing, but 
are not as robust as those from E and CE on GN20. 



AzTECl is not resolved by the SMA in EXT configura- 
tion: its visibility function is flat out to ~ 250 kA, which 
suggests a characteristic angular scale of ^0.5 arcsec 
(see Figure ISj. The inferred flux density from the EXT 
track (-FggoAtm = 13.8 ± 2.3) is furthermore consistent 
with that frorti the C OM track (Fggo^m = 15.6±1.1 mJy; 
lYoimger et al.ll2007f ). The EXT detection is coincident 
with its compact i-band coun terpart in ACS imag ing, a 
faint IRAC 3.6-8 /nm source (j Younger et al.ll2007f l. and 
roughly so with the radio counterpart (-F20cm = 40 ± 13 
/xjy; '^ 0.4 arcsec away) to within the uncertainties - 
though Figure [2] appears to show a potentially signifi- 
cant offset between the radio and SMA positions, this 
is roughly within the total uncertainty in the measure- 
ment of their relative position^. It is not detected i n 
the deep COSMOS 24 ^m imaging (j Younger et al.ll2007f ). 
The submillimeter size of AzTECl is also consistent 
with its 20cm counterpart, which is compact compared 
to the 1.5 X 1.4 arcsec VLA beam. As with GN20, 
the optical counterpart is a B-band dropout, which 
suggests a redshift range consistent with the observed 
radi o-to-submillimeter a nd 24/im-to-subniillimeter ra- 
tios (jYounger et al.ll2007[ ). 

The visibility function for AzTECl and J 1008+063 de- 
rived from the VEX track are shown in Figure |4l Some 
decorrelation on longer baselines - likely the result of 
residual baseline errors in combination with atmospheric 
effects - results in artificial structure in the visibility 
function of J1008-t-063. A gaussian fit to this visibility 
data yields a source size of (0.09±0.02) x (0.07±0.02) arc- 
sec, which describes the effective seeing size for the track, 
and thus the minimum source size which is meaningfully 
probed by these observations. The visibility function 
for AzTECl shows some marginal evidence of being re- 

''' The total uncertainty in the offset between the SMA and VLA 
positions is the combined error from the statistical and systematic 
uncertainties in the SMA position (ctsma stat ~ 0.08 arcsec and 
fSMA.sys ~ 0.05 arcsec), the systematic uncertainty in the test 
quasar's absolute position (cvLBA.sys ^ 15 mas), and the statisti- 
cal and systematic uncertainties in the radio posi tion (cyLA.stat ~ 
HWHM/(S/N) ~ 0.25 arcsec; crvLA.sys $ 55 mas: ISchinnerer et all 
120071) . "This yields a total uncertainty of atot ~ 0.4 — 0.5 arcsec 
which is comparable to the observed discrepancy in Figure |2] 



solved on scales significantly larger than this lower limit: 
a Gaussian fit to this visibility data yields a total flux of 
16.0 ± 5.0 mJy - consistent with the COM, EXT, and 
COM+EXT fits - with a size of (0.30 ± 0.15) x (0.20 ± 
0.10) arcsec. Fixing the total flux to the value derived 
from the COM+EXT data marginally improves this size 
measurement to (0.29+0.13) x (0.18+0.10) arcsec. While 
the statistical uncertainty in the position measurement 
for AzTECl also improves to ^ 0.04 arcsec in both a and 
6, because we calibrate using J1008+063 to remove base- 
line errors it is exactly at the phase center and therefore 
does not provide an estimate of the systematic positional 
uncertainty. Therefore, we quote the position and flux 
derived from the COM+EXT tracks in Tabled 

4. DISCUSSION 

The robust result of these observations is that the far- 
IR emission in both GN20 and AzTECl is small con- 
sidering the very high luminosity of these systems, but 
is clearly extended on ~ kpc scales. This is suggestive 
of mergers as the physical mechanis m driving the bolo- 
metr i c luminosity of these systems (jMihos fc Hernquisfl 
119941 : iHopkins et al.l 1200 6). However in general, and in 
particular for the case of AzTECl, this does not require 
that the far-IR luminosity is contributed only by the 
starburst. Indeed, a significant fraction of the far-IR 
could arise from a dusty torus associated with an ac- 
tive nucleus , which i s generated on significantly smaller 
scales (e.g.. lUrrv fc Padovani 1995). Ideally one would 
like significantly improved resolution continuum imaging 
and resolved gas kinematics via molecular spectroscopy 
to constrain the structure of these sources and dynamical 
state of the star-forming gas in detail - measurements 
which are beyond the capabilities of current facilities 
(e.g., SMA, CARMA, PdBI) but in the near term future 
will likely be accomplished with relative ease by ALMA. 
Nevertheless, if we assume that GN20 and AzTECl are 
sta rburst dominated - as t he typ ical SMG is thought to 
be (| Alexander et al.l 120051 I2008D - and make a series of 
admittedly crude but arguably reasonable assumptions 
about their morphology and kinematics, we find a pre- 
liminary indication that they may be radiating close to 
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Fig. 3. — Top: The real visibility amplitudes as a function oi u — v distance - which shows flux density as a function of decreasing 

angular scale - for GN20 (left) and AzTECl (right), combining all available data from COM and EXT tracks. The real part of the u v 

dat a was binned and scalar averaged, with error bars to indicate the dispersion in a given bin. Also shown are the median radio sizes from 
ICh apman et al. (red, C04; 2004) and Biggs & Ivison (blue, BIOS; 2008). GN20 shows evidence of being partially resolved on scales of ^ 200 
kA - consistent with the median radio sizes of BIOS — while AzTECl is unresolved out to ~ 250 kA - from which we infer a maximum 
angular scale of < 0.25 arsec. The soUd Une indicates a point-source fit to the visibilities (COM+EXT for AzTECl, just COM for GN20), 
with the uncertainty indicated by dashed lines. BOTTOM: The real visibility amplitudes as a function oi u — v distance for test quasars 
observed during EXT tracks for GN20 (J1302+578; right) and AzTECl (JlOOS+063; right). This serves as an empirical verification of 
the phase transfer and source structure. That both test quasars are unresolved by the SMA, with fiat visibility functions, rules out phase 
errors or seeing effects as artificially imposing structure on our targets. 

diameter - in kpc, (T400 is the line-of-sight gas velocity 
dispersion in units of 400 km s ~^, and kioq is the dust 
opacity in units of 100 cm^ g~^ (jMurrav et al.ll200"5l ) . As 
an upper hmit on SFRmax, we adopt kiqo ~ 1, but note 
that many dust models allow for orders of magnitude 
higher opacity, particularly for the ultraviolet radiation 
produced by youn g massive stars during a starburst (e.g., 
iLi fc Draindl200lh . 

Our high resolution observations of AzTECl and 
GN20, which provide a measurement of both the far- 
IR luminosity and physical scale of the starburst re- 
gion in two of the most luminous systems known, of- 
fer an ideal test-bed for this hypothesis. Both targets 
are B-band dropouts (see Figure [2] and §[3]), which in- 
dicates that they likely lie within a relatively narrow 
redsh i ft range 3.5 ^ z ^ 4 .5 (see e.g., ISteidel et al.l 
Il999t iGiavalisco et all l2004bD . While these colors are 
possibly consistent with a very dusty z '^ 2 source, 
the lack of a bright MIPS 24 /im counterpart - aris- 
ing from redshifted ^ 8/im PAH emission - and their 
observed radio-to-submillimetcr flux density ratios - as- 
suming a far-IR SED similar to Arp 220 - are consistent 
with this higher redshift interpretation (jYounger et al.l 
|2007( ). Furthermore template fitting to both the opti- 
cal and far-IR SEDs of AzTECl independently yield 
a consistent result of z ~ 4 (Yun ct al., in prepara- 



or at the Eddington limit of their starburst. 

It has been suggested that feedback from ongo- 
ing star formation sets a physical lim it on the mini- 
mum size of a star forming re gion ( E lmegreenl Il999l : 
iMurrav et all 120051 : iThompson e t al. 20(3). Owing to 
the significant opacity of dust to the ultraviolet light 
produced by young stars, radiation pressure from high 
luminosity star formation regio ns can produce strong 
momentum - driven winds (e.g., iNetzer fc Elitzun Il993t 
lElitzur fcl vczic 2001). These winds are confined by the 
gravitational potential, which for an isothermal sphere 
scales as $ ~ /gCr^logD - where fg is the gas fraction, 
(7 is the stellar velocity dispersion, and D is the diam- 
eter of the starburst region. In the optically thin limit 
- which is appropriate for optically thick clouds with a 
small volume filling factor embed ded in the diffus e inter- 
stellar medium - and assuming a ISalpeteii (|1955f ) initial 
mass function (IMF^), this leads to a maximum star for- 
mation rate (SFR) of: 



SFRn 



900 o-400-Dkpc'- 



■100 



M, 



© yr 



(1) 



where Dkpc is the characteristic physical scale of the star- 
burst - measured via either the Gaussian FWHM or disk 



* Using a IKroupal 1I2OOII) or [Chabrierf 112003) IMF will tend to 
lower SFRmax by ~ 40% l|KennicutH 119981 : IBelll l2003l : IBell et al.l 
[20051) . 
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Fig. 4. — The real visibility amplitudes as a function oi u — v distance for AzTECl (left) and the test quasar J1008+063 (right) using 

data from the VEX track. The real part of the u v data was binned and scalar averaged, with error bars to indicate the dispersion 

in a given bin. J1008+063 shows some evidence of partial decorrelation on the the longer baselines, likely a result of residual baseline 
errors combined with atmospheric effects (see §[2] for a description of the calibration of this track). The results of a point-source fit are 
indicated by the solid line, but its visibility function is well-described by a Gaussian of size ~ 0.08 arcsec (dot-dot-dashed line), which sets 
a lower limit on any meaningful source size derived from these data due to this effective seeing. AzTECl shows evidence of being resolved 
at ~ 500 kA, or angular scales of ^ 0.2 arcsec (dot— dot— dashed line); significantly larger than this seeing size scale and thus a meaningful 
measurement of the source structure. Its total flux derived from just this track is consistent with the point— source fit from the COM+EXT 
tracks (solid line, with dashed lines indicating the uncertainty). 

including tT4oo- Unfortunately, optical redshifts of the 
requisite precision are not available for these sources. 
Similar observations of other SMGs found typical ve- 
locity dispersions of 0-400 ~ 1 for somewhat lower lu- 
minosity systems (median S'sso/jm ~ 8—11 mJy, or 
LpTR « 2 X 1 0"L,y, ; [Neri et al.l l2gM iGreve et al.l [20051 : 
iTacconi et all 120061 I2008D . This is consistent with the 
observed connection between SMGs and present day 3- 
4 L* early-type galaxies, in combination with the star- 
burst mass fractions in the remnants of gas-rich merg- 
ers ex pected from simula tions and observed in local sys- 
tems (jHopkins et al.ll2008ll . In what follows, we use these 
observational constraints and adopt (T400 = 1 with the 
awareness that this parameter is somewhat uncertain for 
the more extreme systems we are studying. 

GN20 has a very high far-IR luminosity, with 
LF/fl(GN20) « 5 X 1013 Lq and SFR(GN20) « 8000 A/q 
yr~i on characteristic a physical scale of £g(GN20) « 4 
and ^d(GN20) « 7 kpc for a Gaussian and elliptical 
disk model respectively. Adopting (T400 = 1, the cor- 
responding Eddington limits for each source model are 
SFR„iax,G ~ 3600 and SFR,„ax,D ~ 3600 Mq yr-\ This 
very luminous SMG is close to or at the Eddington limit 
for a starburst on those scales. It is also interesting to 
note that this size scale is somewhat extended compared 
to a simple R r^ L^^^ - for a disk geometry - or f? ~ L^^^ 
- for a spherical geo metry - scaling of the starburst size of 
local ULIRGs fe.g.. lDownes fc Solomo^ll998Hlono et all 
I2OO7I see also lono et al. 2008, in preparation). 

AzTECl has a far-IR luminosity 1° of 
LF/i?(AzTECl) « 3 X 1013 Lq and SFR(AzTECl) « 
5000 Mq yr^i on a characteristic physical scales of 
^g(AzTECI) « 1.5 and £d(AzTEC1) « 2.5 kpc for 
a Gaussian and elliptical disk model respectively. 
Adopting (7400 = 1, this is significantly larger than 
the Eddington limit for a starburst on this scale, with 
SFR„,ax,G « 1350 and SFR„iax,D ~ 2250 Mq yr-\ 

^^ This differs fr om the luminosity implied by Fig. 3 of 
lYounger et al.l l|2007l ) because it is derived from the far-IR directly, 
not from the near— infrared. 



Fig. 5. — The angular diameter required for a rest-frame bright- 
ness temperature at cosmological redshift Ti,/(l + z) for both GN20 
(solid line) and AzTECl (dashed line). In the optically thick limit, 
in which the dense molecular clouds have a high volume filling 
fraction, the dust temperature T^j = T;,. For GN20, characteris- 
tic angular scale inferred from a Gaussian and disk model yield 
Tf, fs 6(1 -I- z) and 4(1 + 2;) K respectively, which at 2 ~ 4 suggests 
Ta ^ 30 or T^ si 20. For AzTECl, these models yield T;, ai 8(1 + 2) 
and 6(1 + 2) K, which at z rv. 4 gives T^ » 40 and T^ Ri 30 K. 



tion). For this range, the angular diameter distance is 
roughly constant with redshift, so observed angular sizes 
9 correspond to physical scales oi i ^ 7 (0/arcsec) kpc. 
We can also roughly estimate their total far-IR lumi- 
nosity by assuming an Arp 220 template^, which for 
z ^ 1 gives Lf /H « 2 (Fs90um / mJy) 10^^ Lq (s ee also 
iNeri et a"Lll2003[ ) and - assuming a lSalpeterl (|1955D initial 
mass function (IMF) - a star formation rate of SFR « 
340 (Fsgo^m/mJy) Mq yr''^, owing to the strong nega- 
tive fc-correction in the sub millimcter at high-redshift 
(jBlain fc Longaiiill993t iBlain et al. 2002 ). 

CO spectroscopy offers the best route to a measure- 
ment of the dynamical state of the star forming gas, 

^ The inferred luminosity is known to be uncertain by a factor 
of ~ 2 — 3 due to variations in the dust temperature and emis- 
sivity (e.g.. [Blain et aLl 120031) . For example, if we assume a Mrk 
231 template with significant AGN contribution to the far-IR, the 
inferred far-IR luminosity and SFR wi ll be ^ 40 — 60% lower (e.g., 
IStevens et "alll2005l : IHuang et ani2007l ) . 



Two- 


TABLE 3 
-Component Fitting Results 


Name 


-'^1,890/jm -^^2,890^111 ^^ 
[mjy] [niJy] [arcscc] 


GN20'' 
AzTECl^ 


13.5 ±2.4 9.6 ±2.4 0.6 ±0.13 
8.9 ±2.2 4.6 ±2.1 0.3 ±0.11 



'^ Flux of the first component derived from fit- 
ting a two-component point-source model to the 
calibration visibilities.^ Flux of the second compo- 
nent derived from fitting a two-component point- 
source model to the calibration visibilities.'^ The 
separation of the two fitted components, in- 
cluding both the statistical and systematic po- 
sitional uncertainty (see § |2] for details (see also 
lYounger et al.l lSbdT).'^ Derived from a fit to the 
COM+EXT data.'^ Derived from a fit to just the 
VEX data. 

Increasing the dust opacity only aggravates the situa- 
tion. However, some SMGs have been observed with 
(J400 ~ 1.5 — 2, which could explain the discrepancy. 
However, under the assumption that the far-IR luminos- 
ity is dominated by a starburst component, even at high 
velocity dispersion AzTECl is close to its Eddington 
limit. 

If the volume filling factor of dense molecular gas is 
close to unity, and therefore the star forming gas is opti- 
cally thick, then - again assuming a[S alpetcr ( 19551 ) IMF 
- the Eddington limit (jMurrav et al.i i2005i) on the star 
formation rate is 



SFR, 



max, thick 



G 



10V9,o.5't|oo Mq yr-1 (2) 



where /g.0.5 is the gas mass fraction in units of 0.5. This 
is an order of magnitude higher than the SFR of GN20 
and AzTECl. However, in the optically thick limit, the 
dust temperature Td and brightness temperature T;, - de- 
fined as lu = Bi,{Tb), where Ii, is the surface brightness 
and Bi, is the Planck function - are equivalent. The im- 
phed angular diameter Q of both GN20 and AzTECl for 
a given brightness temperature at cosmological redshift 
(Ti,/{1 + z)) is shown in Figure [3 For GN20, charac- 
teristic angular scale inferred from a Gaussian and disk 
model yield Ti, « 6(1 -I- z) and 4(1 -I- z) K respectively, 
which at z - 4 suggests Td « 30 or Td « 20. For 
AzTECl, these models yield Tf, « 8(1 + z) and 6(1 + z) 
K, which at z ~ 4 gives Td « 40 and T^ « 30 K. These 
are all somewhat lower than would be e xpected from 
the temperature-luminosity relation a t low (i Punne et al.l 
I2OOOI: iKlaas et al.ll2001l: I Yang fc Phillips 2007 ), interme- 
diate ("Yang et al."2007^, and high redshift (Blain et al. 
[2003; Chapman et al. 2005; Kovacs et al. 2006). How- 
ever, the brightness temperature represents a lower limit, 
as the inferred dust temperature will increase as the 
opacity r^ decreases - i.e., Td = Ti,{l + z)/(l — e^'^") 
- or if the volume filling factor of optically thick clouds 
is less than unity. There is evidence that in the c ores of 
local ULIRGs rioo^m ^ 1 for A ^ 100 ^m (|Solomon et al.l 



Il997t) , and that the v olume filling factor o f dense molecu- 
lar gas is « 30-70% (jDownes et al.lll993[ ). By analogy, it 
is plausible that GN20 and AzTECl are intermediate be- 
tween the optically thick and optically thin regimes. Fu- 
ture observations a t shorter w a velen gths - e.g., at 350 fim 
with SHARC-H (|Dowell et al.ii2003D - could constrain Td 
independently, and thus help determine the appropriate 
limit. 

Furthermore, our u ~ v coverage does not exclude a 
multi-component structure for either GN20 or AzTECl, 
in particular one with two compact point sources. The 
results of a fit to the calibrated visibilities for this model 
are summarized in Table [3] These angular offsets cor- 
respond to a physical separation of 4 and 2 kpc for 
GN20 and AzTECl respectively, which are consistent 
with dual nuclear starbursts in a late stage me rger (e.g. 
iMihos fc HernauistlfT99l iHopkins et al.ll2006l) . Should 
higher resolution data from either the SMA or ALMA 
confirm this interpretation, it is possible that each of 
these components is at or close to its Eddington limit 
with CT400 ^ 2. 

5. CONCLUSION 

We present high resolution interferometric submillime- 
tcr imaging of t wo of the brightest - an d therefore likely 
most luminous (jBlain fc Longairl |1993[ ) - high redshift 
starburst galaxies known - GN20 and AzTECl. The vis- 
ibility functions for these sources indicate characteristic 
angular sizes of ~ 0.5 — 1.2 and ~ 0.2 — 0.4 arcsec respec- 
tively. Both are B-band dropout optical sources, which 
indicates a redshift of 3.5 ^ z <, 4.5, and thus these an- 
gular size measurement correspond to physical scales of 
4 — 8 and 1.5 — 3 kpc. Assuming a simple morphology and 
a dynamical state typical of high-redshift SMGs, we find 
preliminary evidence that GN20 and AzTECl are both 
close to the limiting luminosity derived via Eddington 
arguments. If these two sources are indicative of a large 
population of hyper luminous starbursts at high redshift, 
this may have important consequences for models of star 
formation and feedback in extreme environments. 
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